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I. INTRODUCTION 



In electrodynamics, the integral equation method (IEM) is a powerful method for the 
electromagnetic wave propagation jl, 2]. From a microscopic perspective, a bulk material 
can be regarded as a collection of molecules (or atoms), each of which is a scatter and emits 
radiation under the action of an external field. The radiated retarded field and the exciting 
field interact to form the resultant transmitted field in the material. This formulation 
is often expressed as an integral equation and naturally leads to the well known Ewald- 
Oseen extinction theorem which is a basis of the scattering theory The IEM relates 
the macroscopic electromagnetic responses of the material with the microscopic features of 
its constituent units, thus gives much deeper physical insight into the interaction of wave 



with material in contrast to the conventional approach of Maxwell theory 



a a 



At the 



same time, the IEM does not require boundary conditions and is advantageous in some 
situations, especially when it is not sufficient to describe the interaction by the macroscopic 



Maxwell equations [4-7]. This approach 
the tip-sample interaction in near optics 
with semi-infinite 



101 ]. a slab 



ms been widely used to study light scattering |5|], 
, 9 1 and so on. By now, the interactions of waves 

it n 

1121 ] or stratified [13[ isotropic dielectric media have been 
treated by the IEM. 

Recently, one kind of artificial anisotropic dielectric-magnetic materials, metamaterials, 
has attracted considerable attention simulated by its fascinating properties [ijj]. To 
understand the unique characteristics of metamaterials it is necessary and appealing to 
take a microscopic viewpoint. On this issue, the IEM has been successfully implemented to 



predict the electric and magnetic resonances in split-ring resonators 16], s tudy the reflection 



of split-ring resonators 



171 ] and the imaging process of super lens 



19 



18|. and explain the 
201 ] . However, most 



Brewster phenomenon for TE wave associated with metamaterials 
previous work only consider the situation of semi-infinite materials and waves incident from 
vacuum. Then, how to treat the propagation of waves through stratified materials, especially 
for incidence from a material rather than vacuum? The latter problem is important because 
it will result in a completely new form of the extinction theorem, as shown in this paper. 

The purpose of this paper is to use the IEM to investigate the properties of wave prop- 
agation in stratified anisotropic dielectric-magnetic media. Our emphasis is to disclose the 
underlying physics of electromagnetic interaction with matter rather that to obtain the well- 
known phenomenal results. On the basis of the superposition principle, we use Hertz vector 
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formulations of radiated fields to study the microscopic mechanism of propagation. We arrive 
at new derivations of the dispersion relation, Snell's law and reflection/transmission coeffi- 
cients by self-consistent analyses. Moreover, we find two generalized forms of the extinction 
theorem. As an application, we use the IEM to study the wave propagation through a slab 
and discuss the physical process of photon tunneling, which are further verified through nu- 
merical simulations. The results lead to a unified framework of the IEM for the propagation 
of wave incident either from a medium or vacuum in stratified dielectric-magnetic materials. 



II. THE INTEGRAL FORMULATION OF HERTZ VECTORS 



In this section we take the microscopic viewpoint in which the matter consists of molecules 
that react to the exciting field like dipoles and formulate the radiated field in the form of a 
Green's function integral of Hertz vectors. 

Let a monochromatic electromagnetic field of Ej = E^o exp (zkj • r — iut) and H« = 
Hj exp (zkj • r — iut) incident from an anisotropic material filling the semi-infinite space 
z < 0. The x — z plane is the plane of incidence and the schematic diagram is in 
Fig. [TJ Since the material responds linearly, all the fields have the same dependence 
of exp (—iut) which will be omitted subsequently for simplicity. Assume the reflected 
fields are E r = E r0 exp (zk r ■ r) and H r = H r0 exp (ik r ■ r) , and the transmitted fields 
are E t = E t0 exp (ik t ■ r) and H t = H t0 exp (ih t ■ r). The permittivity and permeability 
tensors of the n-th layer are simultaneously diagonal in the principal coordinate system, 

Inside the layers, the incident field drive the constituent unities to oscillate and radiate 
ust as dipoles. The radiated electric fields by electric and magnetic dipoles are decided by 

E rad = V(V • n e ) - ^o/io^ - /i V x ^ (i) 



and the magnetic fields generated by dipoles are 



d 2 U m _ dU e 



Urad = v(v • n m ) - ^v-^r + £ oV X (2) 

Here II e and II m are the Hertz vectors, which can be expressed as integrations of the Green's 
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FIG. 1: Schematic diagram for the reflection and refraction of TE waves in stratified media, 
function, 

n e (r)= / ^G(r - r')dr', (3) 

Jy £ o 

n m (r) = / M(r')G(r - r')dr'. (4) 

Jv> 

P is the dipole moment density of electric dipoles and M is that of magnetic dipoles, which 
are related to the internal fields as P„ = e Xen • E n , M n = Xmn ' H n , where the electric 
susceptibility Xen = (en/ £ o) — 1 and the magnetic susceptibility Xmn = (Mn/^o) — 1- Inside 
the layer, the dipoles produce forward waves as well as backward waves. So, we assume the 
dipole moment densities P and M have the following forms 

P n = P nF exp (ik nF ■ r) + P nB exp (-ik nB ■ r), (5) 
M n = M nF exp (ik nF ■ r) + M B exp (-ik nB ■ r), (6) 



where k nF stands for the forward wave and — \t nB stands for the backward wave III, Il2 |. 
The wave numbers are to be decided by the dispersion relation yet to be derived. The Green 
function is G(r — r') = exp (iki\r — r'|)/(47r|r — r'|). To evaluate the Hertz vectors, we firstly 



represent the Green function in the Fourier form. Then, substituting it into Eqs. fl3]) and (j4]) 
and using the delta function definition and contour integration method 1(], 12], the Hertz 
vectors of the n-th layer can be evaluated as 



_ P n F^exp (ik nF ■ r) ^ exp [iq n F • r + i(k nFz QnZ/^n— l] 
kl F - q\ 2q nz (q nz - k nFz ) 

| exp [ig nB ■ r + i(k nFz + q nz )z n \ \ P nB [ exp {-ik nB • r) 

2qnz(qnz + k nFz ) J £ I ^-B ~ 

| exp [-ab nB ■ r - i(k nBz + b nz )z n _-^ exp [-ib nF ■ r - z(A: to b z - qnz)z n ] \ 
2&nz(^n« "I" k nBz ) 2b nz (b nz k nBz ) J 

Zn-1 < ^ < Z n , (7) 

right to the material 



P nF exp [i(A; nj F 2 - q n z)z n -i] ~ exp [i(/c ni?2 - q nz )z n ] . 
n en = — — — exp (tq nF • r) 

fc ^qnzKHnz ^nFz) 

P nB exp [-i{k nBz + b nz )z n -i] - exp [-i{k nBz + b nz )z n ] 
+ oh 71 i t x ex P Hb nB • r), 

z > z n , (8) 

left to the material 



„ i^nf exp 2t( fc nF2 

"I - qnzjZ n \ exp lif/CnFz "I - Zn—l\ 1 ■ \ 

n e „ = — r exp (tq nB • r) 

^Qnz\Qnz + KnFz) 

P nB exp [-i(k nBz - 6 nz )2; n ] - exp [-i(k nBz - 6 nz )2 n _i] . 
+ ^-77 r v ex P {-iKf ■ r), 

£0 ^0 nz [0 nz — K nBz ) 

Z < Z n -!, (9) 

where we have used the following notations 

Qn-F = k nFx it + q nz i, q n _B = k nFx x. — q nz z, q nz = k — k nFx ] 

Kf = k nBx x + b nz z, b nB = k nBx x - b nz z, b 2 nz = k 2 - k 2 nBx . (10) 

Note that the term with exp (*) equals zero if (*) contains Zj — > ±00, which can be found 
in the calculation. Then the Hertz vectors of arbitrary layers can be found in Eqs. (J?])-©. 
Analogous to H ne , n nm can be obtained by interchanging P nF /eo and P nB /eo with M n ^ 
and M nB , respectively. Inserting H ne and U nm into Eqs. (pQ) and the radiated fields can 
be calculated out. 
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III. RADIATED FIELDS AND EWALD-OSEEN EXTINCTION THEOREM 



In this section, we deduce the expressions of radiated fields generated by constituent 
unities in the propagation of waves in stratified anisotropic dielectric-magnetic media. Then 
we use the superposition principle to show how all the fields are combined to yield the 
reflected and transmitted fields inside the individual layers. By a self-consistent analysis we 
obtain the dispersion relation, Snell's law and Fresnel's coefficients. Moreover, we find two 
generalized forms of the extinction theorem. 



A. Fields in the n-th layer 

Inside the n-th layer, the layer by itself and the other layers all radiate under the action 
of external fields. Their radiated fields add up to the real fields inside this layer, that is, 

Enfe/'.O I T?\left . -pn , T7\right . -^right.N+1 /-. -, \ 

n ~ &rad + ^rad + ^rad + &rad + ^rad ■ \ Li ) 

Here the radiated field from the left semi-infinite space 

E *e/t.o _ 1 E *o ex P (* k i ' r )> if z < is vacuum; 

rad ~~ I Q(qoF.-PoF) _ Q(-b 0B ,P 0B ) Qr elge ^ > 

the radiated field generated by the layer itself 

(k n Fz Qnz) z n -l] ^ , D \ 

^rad ~ -p _ 2 + O n ( n _ U I S , f nF ) 

^nFz 'inz ^Hnz\Hnz ^nFz) 

exp[i(k nFz + q nz )z n ] exp [-i(k nBz + b n )z n -i] . 

H ~ 7 —7 x <4iqnB, ir nF ) H — — — — r ^4{-O nB , fnB) 

■^ n nz\ n nz ~T KnFz) ^0 nz [0 nz + KnBz) 

. Q(k„B, Pub) . exp[-i(k nBz -b nz )z n ] 

H 777 77 i \ H{~"nF,-t^nB), U<jJ 



^Bz 2b nz (b nz k n B 

the radiated field generated by left layers (1, • • • ,n— 1) 

rfeft V" ex P WjFz ~ qjz)zj-i} - exp [i(k jFz - q jz )zj} 

| ^ exp H^-fe + 6j Z )^--i] - exp [-ijkjBz + M^j] ^ b ^ p ^ ^ 
■ =1 2bj z (ki z + kjBz) 

the field radiated by right layers (n + 1, • • • , iV) 

™m v- exp [t(fc jFz + q jz )zj] - exp [t(fc jFz + q jz )z j _ l ] 



V- exp [-t(fcj Bz - - exp [-t(fc 3 - Bz - b jz ) Zj ^] wi ^ 

+ 26,(6,, - *,„,) Q(-V,^b). (15) 
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the radiated field from the right semi-infinite space 

v right.N+l -exp [l{k {N+ i) Fz + q{N+l)z)zN) ( . , . 

^rad 71 ^Mlq(AM-l)B,-HiV+l)Fj, UOJ 

^9(JV+l)zi£(AT+l)Fz + <7(JV+l)zJ 

and Q is an auxiliary function 

Q(K, P) = — [K x K x P + (K 2 - £ /i w 2 )P + w/io^oK x Ml exp (iK ■ r) (17) 

where M = Xm 4 [k x (P/e Xe)/ lA/w. There are two cases for the layer. 

Case 1. The layer is an anisotropic dielectric-magnetic slab. Using E n = 
[P nF exp (zk F ■ r) + P nB exp (— ik B ■ r)]/ (e Xne) and inserting Eqs. (fT2|) - ATS]) into Eq. (fTTI) . 
we come to the following conclusions. 

(i). Checking terms associated with the incident field, we find that q\ F = —his = qoF, 
that is, qi z = b\ z and ki Fx = —k\ Bx = ki X that is just the Snell's law. The incident field is 
related to the radiated field as 

exp [i(k {n _ 1)Fz - q( n -i)z)z n -i] 



2<7(n-l)z(<?(n-l)2 — k(n-l)Fz) 

exp [-i{k (n _ l)Bz + q( n -i)z)z n -i\ 



n— \ )z ~r <*(n— \)Bz) 



Q( < l(n-1)F) jP(n-l)F) 
Q(q(n-1)F, P(n-l)s) 



exp [z(A; nF2 - g n2 )2 n _i] exp H(£ nB * + Qw)zn-i] ~, „ x n 

■^4(,qnF, *nF ) H ^ 7 ~ J n ^(qnF, -TnBj = U 



2<7nz('?nz k nFz ) ^qnziqnz ~\~ k nF jz 

n = l,---,N. (18) 

This equation is actually the forward expression of the Ewald-Oseen extinction theorem. It 
shows quantitatively that the forward part of the radiation fields at a location in the layer 
extinguishes those from the regions left to that location. 

(ii). Examining terms of the phase factors exp (iq n 2 • r) and exp (— ih n \ ■ r) in Eq. (TTTI) . 
we know that q„£ = —h nF . At the same time, we have 

exp [i{q( n +i) z + k {n+ i )Fz )z n )} 



2q(n+l)z(q(n+l)z + k 

(n+l)Fz 

exp [-i(k( n+ i) Bz - q(n+l)z)Zn 



2<Z(n+l)z(<?(n+l)z — k(n + \)Bz 



Q(q(n+1)B, P(n+1)F) 
Q(q(n+1)B, P{n+l)B) 



exp[i(q nz + k nFz )z n ] exp [-i(k nBz ~ D v n 

-^ilqn-B) -^n-Fj H ~ 7 J s HKQnB, -t^nB) — U, 



^qnz{qnz ~\~ k nFz ) ^qnz{qnz k nF jz) 

n = l,---,N. (19) 

Eq. ([19]) is the backword form of the Ewald-Oseen extinction theorem. It indicates the 
extinction of the backward vacuum fields in the layer by those from right to the layer. 



Eqs. ([TBI) and (|T9|) are two new forms of the extinction theorem that generalize the previous 



13| 



results to the case of incidence from a material rather than vacuum |4, 

(hi). From terms with the phase factors exp (ik nF ■ r) or exp (— ik nB • r) in Eq. (TTT|) 
follows the dispersion relation 

P P 



for TE waves. So /c„i?2 = k n B Z - Eq. f )20|) relates all the propagation modes in different layers. 

Case 2. The layer is vacuum. Assuming the field in the layer E n = E„i?exp {ik. n F • r) + 
E n B exp (— ik n B ■ r), we arrive at 

exp [i(k(n-i)Fz - Q(n-l)z)Zn-l\ , D v 

&nF — ; : r-\4(q(n-i)F,"(n-l)F) 

■*q(n-l)z{q(n-l)z — K(n-l)Fz) 

exp [-i{k {n -X )B z + q(n-l)z)Zn-l] , , , , 

M\FL{n-l)F,P(n-l)B), (21 J 



2g(n-l)z(0 , (n-l)2 + k(n-l)Bz) 



and 



exp [z(g (n+1)2 + A; (n+1)F2 )2;„,] 

— — - ( — ; r-K4{q(n + i)B,Jr(n+l)F) 

^<?(n+l)z W(n+l)z + K (n+l)Fz) 

exp [-i{k {n+1)B z - q(n + l)z)z n ] , , , , 

M\Si{n+X)B,^{n+l)B) (22 J 



29(n+l)«(9(n+l)* — k( n+ i)Bz 

From Eqs. (121]) and (122]) . it seems that the field in the layer was only related with the 
adjacent regions. However, it is only a mathematical impression because physically all the 
fields from the whole space contribute to the field in the layer, as shown in Eq. ( TTTT) obviously. 

B. The Reflected Field 

In region z < 0, the contributions of radiations from all regions form the real fields, 
yielding 

Ej + E r = Ka d + + Ka 9 d ht - N+1 - (23) 

The field radiated by the material in the region z < can be obtained like Eq. f|T3|) 

go _ Q(koF, Pqf) Q(qoB ; Pqf) QQtpB, Pqb) + Q(— bpF, Pqb) 

^ofz ~ viz 2Qoz(qoz + koFz) k\ Bz — h\ z 2b 0z (b 0z — k 0Bz ) 

Then we have the extinction of backward vacuum waves 

Q(q OB ,PoF) Q(qoB,-PoB) Q(cLib,Pif) Q(qis,Pis) „ 

+ TTTTZ T, \ - o„ I T, V " TZ T, V = °- ( 25 ) 



2qoz(qoz + k Fz) 2qoz(qoz — koB Z ) 2g lz (g lz + k 1Fz ) 2qiz(qiz — ki B z) 
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The incident and reflected fields have the following forms 

P P 

Ej = exp (ik 0F ■ r), E r = exp (— ik Q s ■ r), z < 0. (26) 

By solving the above equations one can obtain the dispersion relation like Eq. ( 120]) and the 
reflected field. If the region z < is vacuum, then 

E r = Q(**,y ■ Q^Pib) z<0 (27) 

2fc(fc + k lFz ) 2q lz (q lz - k lBz ) 



Eq. (|27p is well known in previous work in the case of incidence from vacuum 12|, [13]. 
Differently, Eq. (125|) is a completely new form of the backward extinction theorem because 
the region z < is not vacuum but a material in which the scatters are induced to generate 
radiated field. Using it one can treat the problem of electromagnetic interaction for the 
incidence from either vacuum or a material. 



C. The Transmitted Field 

In the region z > z^, the radiated fields generated by all regions form the transmitted 
field. That is, 

v t = KT + K1t + K + d \ (28) 

where the radiated field produced by dipoles in the region z > zn is 
v n+i Q(<1(n+i)f,P(n+i)f) , exp [i{k {N+1)Fz - q {N+1)z )z N ] 

Kad - -p -3 7z 1 rviiq(jv+i)F,-njv+i)Fj- i^yj 



Hence, one can easily obtain 

exp [i(k (N+1)Fz - g( N+ i) z )z N ] 

2<l{N+l)z(<l(N+l)z ~ k{N+l)Fz) 



Q(q(7V+l)F, P(N+1)f) 



exp [i(q Nz - k NFz )z N ] exp [-i{k NBz + q Nz )z N ] 

"WVQjvb, "jvfJ o 7 — ; s HKSInb, ^nb) = U 



2Qnz(<1nz ~ k NFz ) 2qNz(<lNz + k NBz ) 

(30) 

which describes that the forward vacuum wave generated by the material in the right self- 
half space is extinguished by the counterpart from to the region z < Zn- The transmitted 
field is assumed to be 

E t = — (JV+1)F exp (ik {N+ i )F ■ r) (31) 

£0X(N+l)e 
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If the region z > zn is vacuum, then 



exp[i(k NFz -q Nz )z N ] exp[-i{kNBz + qNz)z N ] r . f „ s 

— ^ / -j \ ^(.qjVF; *^NF) 7, 7 n \ HK^iNF-, -t^NB), 

2qNz{qNz ~ k NFz ) 2q Nz (q Nz + k NB z) 

z>0. (32) 

By solving the equations one can obtain the transmitted field and the dispersion relation in 
the region. 

D. Transformation Matrix 

Actually, it is not difficult to obtain the reflected field and the transmitted field. Using 
Eqs. (IISI). pgjl , (l25]h and pi (Eqs. (|2T]), d22]) , (1271). and (132} if the associated region is 
vacuum), one can derive the mutual relationship of E n in all regions 

= L n{n+1) I E[n+1)F I , n = 0, • • • , {N + 1). (33) 

y ^(n+l)B / 

The transformation matrix L n ( n+1 ) for TE waves is 




L n (n+1) 



exp [i(A;(„+i)F 2 - k nFz )z n ] R n (n+i) exp [-i(fc njF2 + k {n+1)Fz 
T n{n+i) \ R n(n+l) exp [i{k nFz + fc (n+ i) Fz )z n ] exp [i(k nFz - k {n+1)Fz )z n ] 



(34) 

where 

D ^{n+l)xk n Fz ^nxk(n+l)Fz rp ^^{n+^xknF z /oe:\ 

-Kn(n+i) = -j : ; , J-n{n+i) = ; : ; ■ (3oJ 

H{n+l)xKnFz T ^nx^(n+l)Fz f^(n+l)x^nFz T ^nx^(n+l)Fz 

Note that E(n+i)b = 0. Then, with Eqs. (1261) . (13T1) and ( 1331) we finally obtain the real 
reflected and transmitted fields in individual layers. 

The magnetic fields can be obtained just as the corresponding electric fields after replacing 
the function Q(K,P) with 0(K,M), 

0(K, M) = -(KxKxM + K 2 M - e {i u 2 M -wKxP) exp (iK ■ r), (36) 

where P = e Xe • E = -e Xe ■ [k x (M/x m )/e]/u;. 

By now, we have derived in a new way the reflected/transmitted fields and the reflec- 
tion/transmission coefficients in stratified media. Although involving some seemingly com- 
plicated integrals that nevertheless have simple solutions, the IEM leads almost miraculously 
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to the correct results. Obviously, Eqs. (|33|) and (|37|) are in agreement with the results ob- 
tained by the formal approach of Maxwell theory |l|. 

It seems from Eqs. (ITS]) . (]T9]) . or (I5U|) as if the extinction only occurs on the interfaces 
4|. In fact, however, the vacuum fields are extinct at every points in all regions, as shown 
evidently in Eqs. (TIT]) . (]2U]) or (T5T|) . In the above analyses of this section we do not presume 



any boundary condition. Instead, one can 



boundary conditions in Maxwell theory 



ind that the result of Eq. (]33j) plays the role of 



23]. 



IV. PROPAGATION THROUGH A SLAB 

As an application of the IEM, we discuss the case of a three-layer configuration. For the 
case of N = 1, by Eqs. (126]) and f]3T]) we obtain the reflection coefficient Rte = E r0 /E iQ and 
the transmission coefficient Tte = E t0 /E i0 for TE waves 

= gpi + R 12 e l2k ^ _ T 01 T 12 e ik ^ 

TE 1 + R 01 R 12 e i2k ^ ' TE 1 + R Q1 R 12 e i2k ^ ' 1 ' 

To further demonstrate the validity of the IEM, we calculate the propagating fields in photon 
tunneling and compare the results by IEM with those by solving Maxwell equations plus 
boundary conditions Q]. 

Consider a plane wave incident through a slab sandwiched between two prisms at the 
angle larger than critical angle. We calculate the transmitted field in the second prism by 
the IEM and the Maxwell approach, denoted by |E t |j#M and \E t \ Maxwe u respectively, as 
shown in Fig. |2] From that figure, we can see the two results overlap each other, indicating 
that the IEM is consistent with Maxwell approach and thus is valid. In order to disclose the 
physical mechanism of photon tunneling, we calculate the radiated fields |E^ rce | from the 
left to the second prism, which is the source of the transmitted field. One can see that the 
amplitude of E^ rce decreases exponentially as the width of the slab increases. Such a trend 
also occurs to \E t \ IEM . This phenomenon can be understood easily as the transmitted field 
results from the radiated field induced by E^ rce . Especially, we find that when |E^ rce | 
approaches zero, the transmitted field is also close to zero. In other words, photon tunneling 
can not occur without the radiated field E^ rce from the left region. Actually, detailed 
analyses show that the physical process of photon tunneling is that the vacuum field E^ rce 
induce the dipoles in the second prism to radiate two kinds of waves; one of the radiated 
fields extinct E^ rce and the other forms the final transmitted field, leading to photon 
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FIG. 2: The field magnitude as a function of the width of the slab d. A TE-polarized plane wave 
is incident on the air gap sandwiched by two prisms with the dielectric constant 2.576 at the angle 
vr/4. 

tunneling. Based on the above results, we make numerical simulations of photon tunneling 
for slabs with different widths in Fig. |3j For completeness, we give results for both positive 
and negative indices. From Fig. |3l one can see that the photon tunneling mostly occurs as 
the slab is thinner than a wavelength and the transmitted field increases as the slab width 
decreases. At the same time, positive or negative Goos-Hanchen shifts are confirmed in 
photon tunneling as observed in j^]. 



V. CONCLUSIONS 



In summary, we have extended the IEM to investigate the propagation of electromag- 
netic waves through stratified anisotropic dielectric-magnetic media. From a microscopic 
viewpoint in which the matter consists of molecular scatters, we used the superposition 
principle and the integral formulation of Hertz vectors to analyze the interaction of field 
with matter. We presented new derivations of the dispersion relation, Snell's law and the 
reflection/transmission coefficients by self-consistent analyses. Applying the IEM, we in- 
vestigate the wave propagation through a slab and discuss the physical process of photon 
tunneling. Verified by numerical simulations, the results are in agreement with those ob- 
tained by macroscopic Maxwell theory and disclose the underlying physics of wave propaga- 
tion through stratified materials, which is the emphasis of the present paper. Moreover, we 
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FIG. 3: Field distributions for a Gaussian beam through slabs with different widths d. (a) d = 2A, 
(b) d = A, (c) d = 0.5A, (d) d = 0.1A for the slab with e± y = 1 and fj,± x = fii z = 1; (e) d = 2A, (f) 
d = A, (g) d = 0.5A; (h) d = 0.1A for the slab with e± y = —1 and fii x = fj,± z = —1. The slab is 
sandwiched by two prisms with the dielectric constant 2.576 and the incident angle is it /A. 



found two new forms of the extinction theorem generalized to incidence not from vacuum 
but from a material, which may enrich the theory of molecular optics We also 

demonstrated that the_extinction occurs at an arbitrary location inside the materials, not 
just on the interfaces 



refraction index 
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4j. Obviously, we recovered previous results for materials of positive 
More importantly, our results lead to a unified framework of the 
IEM for wave propagation in stratified dielectric-magnetic materials. 

The unified IEM can be applied to other problems of propagation through dielectric- 
magnetic materials, including metamaterials. This is because the units of metamaterials, 
whose dimensions are much less than the wavelength, react to exciting fields just as molecular 
scatters in conventional media [22]. Especially, the method can treat propagations of wave 
incident either from a medium or vacuum, in contrast with previous work [3, [l^] in which 
only the incidence from vacuum is considered. Besides, the same framework not only enables 
the handling of arbitrary shaped objects made of continuous matter, but allows one to treat 
discrete particles. The associated coupling can be included in the IEM without any formal 
difficulty j^l, 9]. It is hoped that the framework and the microscopic perspective of the IEM 
may enable more unique properties in metamaterials to be understood and more applications 
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to be developed. 
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